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We perform large-scale Monte Carlo simulations using the Machta-Newman-Chayes algorithms 
to study the critical behavior of both the diluted antiferromagnet in a field with 30% dilution and 
the random-field Ising model with Gaussian random fields for different field strengths. Analytical 
calculations by Cardy [Phys. Rev. B 29, 505 (1984)] predict that both models map onto each other 
and share the same universality class in the limit of vanishing fields. However, a detailed finite-size 
scaling analysis of both the Binder cumulant and the two-point finite-size correlation length suggests 
that even in the limit of small fields, where the mapping is expected to work, both models are not in 
the same universality class. Therefore, care should be taken when interpreting (experimental) data 
for diluted antiferromagnets in a field using the random-field Ising model. Based on our numerical 
data, we present analytical expressions for the phase boundaries of both models. 
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I. INTRODUCTION 

The random-field Ising models (RFIM) is of 
paramount importance in the field of disordered 
systemsi^"— A plethora of problems across disciplines can 
be studied via the RFIM, ranging from the thermody- 
namics of disordered magnetS)^ hysteresis in magnetic 
systems and Barkhausen noise^— tunable domain-wall 
pinning,— the random pinning of polymers,— and even 
water seepage in porous media. As such, the RFIM is 
still under intense theoretical, as well as numerical and 
experimental scrutiny. 

More recently, the RFIM has been realized in di- 
luted dipolar magnets in a transverse field such as 
LiHoxYi_xF4. However, most experimental studies fo- 
cus on diluted antiferromagnets in a field (DAFF), such 
as FexZni_xF2i^'i^"— Fishman and Aharonyi^ were the 
first to note that the DAFF can be described by the 
RFIM and Cardyii predicted that the critical behav- 
ior of both models should be in the same universality 
class in the limit of small fields. The work of Fish- 
man and Aharony^ as well as Cardyii therefore opened 
the door for intense experimental investigation of the 
RFIM via DAFF materials. However, early experiments 
and simulations already hinted towards discrepancies be- 
tween experimental and numerical estimates of the criti- 
cal exponentsi^ However, exact ground-state calculations 
using moderate system sizes suggested an agreement be- 
tween the critical exponents for both models.— 

In this work we perform detailed Monte Carlo sim- 
ulations of both the RFIM and the DAFF. The latter 
is studied at 30% dilution, i.e., below the percolation 
threshold. Using a finite-size scaling analysis of both the 
Binder cumulant and the two-point finite-size correlation 
function, we show that even in the limit of small fields — 
where the Cardy mappingi^ is expected to work — both 
models seem to be in different universality classes. There- 



fore, care should be taken when making predictions for 
the critical behavior of the RFIM using experiments on 
DAFF materials. Finally, we present heuristic analytical 
expressions based on our numerical data for the phase 
boundaries of both models to help guide experimental 
studies. 

The manuscript is structured as follows. In Sec. [Til we 
introduce both the RFIM and the DAFF, followed by an 
explanation of the used algorithms in Sec. Illli as well as 
the measured quantities in Sec. IIVI In Sec. |V] we show 
our numerical results, followed by a detailed discussion 
of the phase boundaries and universality between both 
models in Sec. I VII 



II. MODELS 

The Hamiltonian of the diluted antiferromagnet in a 
field (DAFF) is given by 



-\-J EiEjSiSj — B EiSi , 



(1) 



and the Hamiltonian for the random-field Ising model 
(RFIM) is 



T~Lt, 



-J ^ SiSj 



(2) 



In Eqs. ([T|) and ^ Si G {±1} represent Ising spins, J = 1 
is the coupling constant between two adjacent spins, and 
denotes a sum over nearest neighbors. The linear 
term in Si couples to an external field: For the DAFF 
it is an externally-applied uniform field B, whereas for 
the RFIM the spins couple to a random field of strength 
hSi, where the Si are quenched random variables chosen 
from a Gaussian distribution with zero mean and stan- 
dard deviation unity. This means that the typical field 



2 



has strength h. In the DAFF Si G {0, 1} represents the 
site dilution, where each site is randomly and indepen- 
dently occupied by a spin (e^ = 1) with probability p. 
Here, we fix the dilution to 1 — p = 0.3. Both models 
are studied in three space dimensions on a lattice with 
N = spins, L the linear size of the lattice. 



III. ALGORITHM 

The simulations are done using the Machta-Newman- 
Chayes replica- exchange (MNC) algorithmic combined 
with single-spin Metropolis Monte Carloi^i^ The MNC 
algorithm is a mixture of the Swendsen-Wang exchange 
algorithm^ and simulated tempering Monte Carloi^i^ 
Note that the latter is not efficient when simulating 
random-field systemsi^ The advantage of the MNC al- 
gorithm over standard parallel tempering lies in the fact 
that we can choose any path in the field-temperature 
plane. Although parallel tempering can also be imple- 
mented with a variable field, the method does not per- 
form efficiently when systems have disorder 

In the MNC algorithmic a cluster of connected spins 
is grown between two replicas with the same disorder 
but at different points in the phase space, i.e., {T,B) 
and {T',B')^ where T represents the temperature and 
B the external field (here for the case of the DAFF). 
Starting from an arbitrary spin with different sign in both 
realizations, adjacent spins pointing in the same direction 
are successively added to the cluster with probability 



p(/3,/3') = l-cxp{-2(/3 + /?')}, 



(3) 



where /3 = l/T is the inverse temperature. Once no 
more spins can be added to the cluster C, it flips with 
the Metropolis probability^! min{l, exp(— E)}, where 

EoAFF = 2sign{C)[{/3-/3'){n++-n^^) + {B-B')\C\\ (4) 
for the DAFF, and for the RFIM 

Ehpim = 2sign(C) [il3 ~ P'){n++ - 7i__) -{h- h') ^ 5, 

iec 

. (5) 

Here \C\ is the number of spins in the cluster, sign(C) 
the orientation of the spin in the replica having inverse 

temperature /3, n++ and n arc the number of bonds 

connecting to nearest neighbors of the cluster with spin 
up and spin down in both replicas, respectively. After 
each cluster update, attempts to flip single spins 

are performed, where L is the linear size of the system. 

As stated before, the MNC algorithm enables us to 
perform simulations along any arbitrary path in phase 
space. We denote such path a rephca chain (RC). The 
phase boundaries for the RFIM and DAFF in the field- 
temperature plane are well described by ellipses (see be- 
low). To reduce corrections to finite-size scaling^S we 
therefore choose paths in the field-temperature plane 



that cut the phase boundaries at an as orthogonal an- 
gle as possible. This means that, in general, T h for 
the RFIM and T ~ B for the DAFF. To ensure efficient 
mixing and therefore fast convergence of the Monte Carlo 
method, we additionally connect the point with the high- 
est field within the disordered phase to another RC that 
runs parallel to the approximated phase boundary to a 
temperature T > and B = (/i = for the RFIM), 
where Tc is the critical temperature of the model at zero 
field (see Fig. [SJ light dashed lines). This end point of 
the second RC is simulated efficiently by the Wolff cluster 
algorithm)^ Simulation parameters are listed in Tables 
U and in for the RFIM and DAFF, for the first RCs, re- 
spectively. 

Finally, we also study the DAFF at zero temperature 
using the method introduced in Ref . Isol. Here, the DAFF 
is mapped onto a grapb^ with N nodes {N is the num- 
ber of spins) attached to a source and a sink node, all 
connected in a distinct manner via edges with positive 
edge weights. The edge weights are calculated depend- 
ing on the local staggered field, i.e., ±B. The maximum 
flow/minimum cut is obtained using the algorithm intro- 
duced in Ref. |32| . The minimum cut is a direct represen- 
tation of the ground-state spin conflguration from which 
derived quantities, such as a zero-temperature Binder ra- 
tio, can be calculated. Note that the method takes the 
ground-state degeneracy into account. The simulation 
parameters for the DAFF at zero temperature are shown 
in Table Uni 



IV. OBSERVABLES 

Both the DAFF and RFIM undergo second-order 
phase transitions as a function of temperature and field. 
To pinpoint the transition temperature we measure the 
Binder cumulant^ as well as the two-point finite-size 
correlation functioui^"— To compute these obscrvables, 
we measure the magnetization per spin 



1 ^ 



(6) 



For the DAFF we measure the staggered magnetization, 
i.e., each second spin is counted opposite to its orienta- 
tion in a three-dimensional checker-board manner. For 
simplicity, we refer to the staggered magnetization also 
as M. An antiferromagnetically-ordcred spin configura- 
tion has therefore M = 1. A Binder cumulant for M can 
then be defined via 



(7) 



where (■ • ■ } represents a thermal average and [■ • ■ ]av an 
average over disorder (field or dilution configurations) for 
a fixed value of h (RFIM) or B (DAFF). Close to criti- 
cality the Binder ratio scales as 



.g(T,i) = G[Li/'^(r-T,)] 



(8) 
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TABLE I: Simulation parameters for the RFIM along dif- 
ferent nontrivial paths of the type h = a + bT in the h-T 
plane for different linear system sizes L (the first two path 
types have 6 = 0). A'^sa is the number of disorder realiza- 
tions. Nt corresponds to the number of temperatures (points) 
along the simulation path. Tmin and Tmax are the lowest and 
highest temperature simulated, respectively. The equilibra- 
tion/measurement times are 2^ Monte Carlo sweeps. 
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where G is a universal function. Note that for the DAFF, 
when T = 0, g{B,L) = C'lL^/^iB - Be)]. To com- 
pute the two-point finite-size correlation function we first 
calculate the wave- vector-dependent susceptibility (along 
the X direction) via 



X(k) = 



TABLE II: Simulation parameters for the DAFF along non- 
trivial paths of the type B = a + bT in the B-T plane for 
different linear system sizes L (the first two path types have 
6 = 0). Asa is the number of disorder realizations. A^r cor- 
responds to the number of temperatures (points) along the 
simulation path. Tmin and Tmax are the smallest and the 
highest temperatures of the RC, respectively. The equilibra- 
tion/measurement times are 2^ Monte Carlo sweeps. 
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The two-point finite-size correlation function is then 
given by 



(9) 



1 



x(o) 



2 sin(fcmi„/2) y 



(10) 



4 



TABLE III: Simulation parameters for the DAFF at zero 
temperature for different fields B and for different linear sys- 
tem sizes L. Nsa. is the number of disorder realizations. Bmin 
and -Bmax are the lowest and highest fields simulated and Nb 
corresponds to the number of fields simulated to perform a 
finite-size scaling analysis. 



B 

max Nb 

24 10302 2.00 4.30 31 

32 2091 2.40 2.70 16 

48 2091 2.10 2.80 17 

64 2091 2.30 2.70 21 

72 2040 2.30 2.54 17 

96 5100 2.30 2.54 17 

128 3586 2.30 2.47 22 



the RFIM are summarized in Table HVl for the DAFF in 
Table m 

TABLE IV: Critical temperature Tc and critical field he com- 
puted from a finite-size scaling analysis of the two-point finite- 
size correlation function for the RFIM. v is the critical expo- 
nent of the correlation length. 
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with kniin = (27r/L, 0, 0). The two-point finite-size corre- 
lation function scales as 

iL{T,L)lL = X{L^I-{T -T,)\. (11) 

Using both the Binder ratio and the two-point finite- 
size correlation function allows us to perform a detailed 
finite-size scaling analysis to determine the critical ex- 
ponent V, as well as to test if both models share the 
same universality class using the method introduced in 
Ref. [S^- To obtain an optimal data collapse we use 
a Levenberg-Marquardt minimization combined with a 
bootstrap analysis, see Ref. 113. This allows us to deter- 
mine the optimal values of the critical parameters and 
V with a statistical error bar by fitting the data to a third- 
order polynomial that approximates the scaling functions 
G[x) and X(x) close to x = 0, where x = L^I''{T - T^). 

V. RESULTS 

The critical parameters for both the RFIM and the 
DAFF have been computed via a finite-size scaling 
analysis of the two-point finite-size correlation function 
[Eq. ([TT]) ] along the different simulation paths. Finite- 
size corrections can be large for small system sizes and 
are strongly field dependent, which is why for some exter- 
nal fields in both models we do not include small systems 
in the finite-size scaling analysis used to determine the 
critical parameters. To illustrate the typical behavior, in 
Fig. [TJ left panel, we show the two-point finite-size cor- 
relation function for the DAFF for B = 1.0 and different 
system sizes. The data cross at a point, therefore sig- 
naling the existence of a phase transition. Note that for 
this particular field corrections to scaling are manageable 
and the data scale well, as can be seen in Fig. [TJ right 
panel. However, this is not always the case, especially 
when the external field is large. For the RFIM correc- 
tions to scaling are considerably stronger, even at small 
fields, see Fig. [H 

Using finite-size scaling we determine the location of 
the critical points, as well as the associated critical ex- 
ponent V for the different simulation paths. Data for 



TABLE V: Critical temperature Tc and critical field Be com- 
puted from a finite-size scaling analysis of the two-point finite- 
size correlation function for the DAFF. v is the critical expo- 
nent of the correlation length. The last line lists data from 
zero-temperature simulations (see text). The estimate of the 
critical field Be is obtained from a finite-size scaling analysis 
of the zero-temperature Binder ratio. 
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To determine the critical field at zero tempera- 
ture for the DAFF we compu te g round states with the 
algorithm introduced in Ref. l30l The same finite-size 
scaling technique as used for the two-point finite-size 
correlation function (see above) can be used to analyze 
the ground-state Binder cumulant. The data collapse is 
shown in Fig. [3l The results for the critical point and 
the correlation-length exponent at zero temperature are 
stated in the last line of Tab. |Vl 

We also determine the peak position of the fluctuations 
of the staggered magnetization of the ground-states: 

T{B) = {[M%, - [M]l) . (12) 

This approach has proven to be quite accurate in pre- 
vious studies for the susceptibility4^ Because the fluc- 
tuations peak at the putative transition, we fit a Gaus- 
sian to the peak and determine its precise location. Er- 
ror bars are determined via a configurational bootstrap 
analysis4i Figure |4l left panel, shows the fiuctuations at 
zero temperature and as a function of the applied field 
B. The peaks are well described by Gaussians. The 
right panel of Fig. 2] shows an extrapolation of the peak 
position to infinite system size assuming the functional 
form Bc{L) = Be + aL~'^ . The best fit is obtained for 
Be = 2.34(2) [bj ^ 1.25(9)], in agreement with the esti- 
mate using the Binder cumulant, see Tab. |Vl 
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FIG. 1: (Color online) Left; Two-point finite-size correlation function £,l/L vs temperature T for the DAFF with B = 1.0 and 
different linear system sizes L. Finite-size corrections are small and the data cross at one point signaling a transition. Right; 
Finite-size scaling collapse of the data in the left panel. The best collapse is obtained with Tc ~ 2.807 and v ~ 1.2. 




FIG. 2; (Color online) Left; Two-point finite-size correlation function £,l/L vs temperature T for the RFIM with h = 0.225 
and different linear system sizes L. Finite-size corrections are large. Right; Finite-size scaling collapse of the data in the left 
panel. Because of the large corrections to scaling, only data for L > 28 are used in the scaling collapse. Data for L < 24 (light 
shaded) are not included in the data collapse and shown to illustrate the corrections to scaling. The best collapse is obtained 
with Tc ^ 4.481 and v ^ 1.39. 



Combining the data in Table IVl with some values from 
the literature^"— we can approximate to high precision 
the phase boundary for the DAFF via 



with f° « 2.980 and B° w 2.31. Similarly, using 
the data from Table IIVI and known values from the 
literature^i^i^ we obtain for the RFIM 
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FIG. 3: (Color online) Data collapse of the zero-temperature 
Binder cumulant of the DAFF as a function of the reduced 
scaling variable {B — Bc)L^^" for different system sizes. The 
best collapse is obtained for ~ 2.32 and v « 1.43. 



with /i2 = 2.27 and T° = 4.5115. Note that the critical 
phase boundary points and hP^ have been determined 
to high precision in the literature. For the DAFF T° and 
are approximated, but agree with the numerical esti- 
mates we present. In Fig. [5] we show the phase boundaries 
for the DAFF (left panel) and the RFIM (right panel), 
together with the simulated critical points. The dashed 
lines represent the simulation paths taken. 



VI. DISCUSSION 

CardjJi predicted an equivalence between the DAFF 
and the RFIM for small applied fields. This equivalence 
is often quoted in experimental studies where materials 
which are diluted antiferromagnets in a field are then 
described using the RFIM (sec, for example, Refs. illl- 
0). 

Equation (15) in Ref. [ij maps the RFIM onto the 
DAFF: 



h{T) 



p(i -p)(TP"^7r)2(B/r)2 



(15) 



Here, p ^ 0.7, T?"™ ^ 4.5115 and 9"''' = 2dJ = 6 is the 
mean-field coupling strength. We can now use the ob- 
tained phase boundaries [Eqs. ([T^ and p^ ] to compare 
both models. Figure [S] shows the phase boundary for the 
RFIM (solid line, the circles represent the obtained crit- 
ical points along the different simulation paths) together 
with the phase boundary for the DAFF mapped onto 
the RFIM space using Eq. ([TS|) (dashed line, the squares 



represent the obtained critical points along the differ- 
ent simulation paths for the DAFF). For random-field 
strengths of up to h ?» 1.2 — which means field strengths 
of up to i? w 1.6 for the DAFF — there is an approx- 
imate correspondence between both models. However, 
as the figure clearly illustrates, strictly speaking the cor- 
respondence only seems to work in the limit of /i — ^ 
{h < 0.3). It is of importance to take these limita- 
tions of the Cardy mappingii into account when study- 
ing diluted antiferromagnets in an external field experi- 
mentally while attempting to describe the data analyt- 
ically using the RFIM. Furthermore, a basic finite-size 
scaling analysis leads to no systematic deviations of the 
correlation-length exponent v. Including the estimates 
for rough simulations at high fields, our results support 



ly = 1.36(14) 



(16) 



for the range of fields studied, in agreement with previ- 
ous studies, such as t'RPiM = 1.37(9)^^ ly = 1.20(5)^^ 
from experiments on Fco.85Zno.15F2 {p = 0.85), or ^ = 
1.40(6) from the disconnected part of the susceptibility of 
Feo.93Zno.07F2 {p = 0.93)4^ Our results are summarized 
in Fig. [T] As can be clearly seen, the difference between 
the estimates for the critical exponent of the correlation 
length for both models is marginal and within error bars. 

However, the error bars are large and therefore a more 
detailed study needs to be performed. To truly discern 
if both models are in the same universality class, in ad- 
dition of having one (apparently) agreeing critical expo- 
nent, one would have to compute a second critical expo- 
nent. Typically, the critical exponent t] that describes 
the decay of the susceptibility at criticality is computed. 
However, and especially for large fields, corrections to 
scaling arc large for both the RFIM and the DAFF. A 
different approach is the computation of the critical expo- 
nent a that describes the divergence of the specific heat. 
However, for both the RFIM and the DAFF a is close 
to zeroj^i^ Therefore, simulations of very large system 
sizes which are currently not accessible numerically are 
required. 

Fortunately, there is a simple yet more sensitive 
method to verify if two different systems share the same 
universality class without having to compute any crit- 
ical exponentsi^2ill Both the Binder cumulant and the 
two-point finite-size correlation function divided by the 
system size are dimensionless quantities. By plotting one 
as a function of the other, nonuniversal quantities cancel 
out<^ For a given system, once large enough system sizes 
are reached such that corrections to scaling are negligi- 
ble, the data for all system sizes collapse onto a universal 
curve within error bars. If two systems share the same 
critical exponent v, we expect that all data should col- 
lapse onto the same universal curve within error bars and, 
in particular, that the estimates of the Binder cumulant 
and the two-point finite-size correlation function agree at 
the putative critical point (s). We therefore would expect 
that data sets of ^(^l/L) for both the DAFF and the 
RFIM should agree for all simulated temperatures and. 
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FIG. 4; (Color online) Left: Fluctuations of the staggered magnetization of the DAFF as a function of applied field B for 
different system sizes. The peak positions signals the presence of a transition. The data are well described by a Gaussian close 
to the peak (solid lines). To determine the thermodynamic critical field Be we extrapolate the data to infinite system size (right 
panel) using Bc{L) — Be + aL~'^ . The best fit is obtained for Be = 2.34(2) and uj = 1.25(9). The red (filled) point represents 
the thermodynamic extrapolation, Be — 2.34(2). 




FIG. 5: (Color online) Left: Empirical phase boundary of the DAFF {p = 0.7). The red point is from Ref. |3^ . the coral point 
from Ref. Isol . and the purple points from Ref. |43. Our data agree within error bars with these previous studies. The blue 
(solid) curve is given by Eq. (|13p . The dashed lines represent the parts of the simulation paths that cross the phase boundary. 
The light shaded line is an example of the second RC that runs parallel to the approximated phase boundary to a temperature 
T > Tc and B = to speed up equilibration. Right: Empirical phase boundary of the RFIM. The zero-field critical temperature 
is Tc — 4.5115^ and = 2.270^ (grey open circles). The red (solid) curve is given by Eq. (|14|) . The dashed lines represent 
the parts of the simulation paths that cross the phase boundary. Again, the light shaded line shows an example of the second 
RC that runs parallel to the approximated phase boundary to a temperature T > Tc and h — to speed up equilibration. 
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FIG. 6: (Color online) Phase boundary of the RFIM (solid 
line, from Fig. [H right panel). The temperature axis has been 
normalized with Tc{h = 0) = 4.5115. The circles represent 
the different estimates of the critical points along the different 
simulation paths. The dashed line is the phase boundary 
computed by applying Eq. (|15|l to the data of the DAFF. 
Squares represent the different critical points simulated for the 
DAFF along the different simulation paths. An approximate 
correspondence between the phase boundaries only works for 
fields h < 1.2 {B < 1.6 for the DAFF). 



in particular for T ^ Tc- 

Figure |8] shows the Binder cumulant function of 
the two-point finite-size correlation function divided by 
the system size for both the DAFF and the RFIM. The 
left set of points (reddish/light tones, circles) are for 
the RFIM. Data for the different simulation paths used 
collapse onto a master curve. The right set of points 
(greenish/dark tones, squares) are for the DAFF. Again, 
all data collapse onto a master curve for all simulation 
paths taken. However, the data sets for the RFIM and 
the DAFF do not agree, except in the trivial limit where 
g{T) 1. The large circles for the RFIM (squares for 
the DAFF) represent our estimates of g{£,L/L) at T = Tc- 
As can be seen, the data for both models do not agree 
(i.e., a large circle should sit on top of a large square), 
something which is even more clear when zooming into 
the boxed area (inset). Note that the large error bars are 
due to the uncertainty of the critical temperature. 



VII. CONCLUSIONS 

We have performed extensive Monte Carlo simulations 
of the diluted antifcrromagnct in a field at 30% dilution 



FIG. 7: (Color online) Critical exponent v as function of the 
field h (RFIM) and B (DAFF). The labels on the upper axis 
correspond to the random-field strength h (RFIM), those on 
the lower axis to the external field B (DAFF). The weighted 
mean is = 1.36(14) (grey line) and the weighted error is 
represented by the shaded (light blue) area. The difference 
between S'daff (blue dashed line) and i'rfim (red dashed line) 
is marginal in comparison to the error-bars of the data points. 
The RFIM ground-state value is taken from Ref. |4^ . 



(p = 0.7) and the random-field Ising model. Using these 
data we show that the phase boundaries for both models 
are well described by ellipsoids (see Fig. [SJ. In addi- 
tion, using zero-temperature heuristic methods, we com- 
pute the zero-temperature critical point for the DAFF 
with 30% dilution {p = 0.7). We expect that the phase 
boundary for other dilutions will be similar, albeit with 
different nonuniversal parameters. 

Furthermore, we numerically study the equivalence of 
the RFIM and the DAFF as predicted by Cardy-il Our 
results show that only in the limit of small fields both 
phase boundaries map onto each other. 

Finally, we perform a finite-size scaling analysis to de- 
termine the critical exponent v of the correlation length. 
Our results from the two-point finite-size correlation 
function suggest that the exponent v agrees within error 
bars for both the RFIM and the DAFF. However, error 
bars are large. To circumvent this problem, we study the 
Binder cumulant as a function of the two-point finite-size 
correlation function divided by the system size and show 
that both models apparently do not share the same uni- 
versality class. Clearly, more detailed simulations need to 
be performed to fully discern the critical behavior of both 
models and fully determine their universality classes. It 
would be interesting to also measure the critical behav- 
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FIG. 8: (Color online) Binder ratio g as a function of the two-point finite-size correlation function divided by the system size 
^l/L for several system sizes and simulation paths. The data for the RFIM and DAFF collapse onto two distinct set of curves, 
suggesting that both models do not share the same universality class. The left set of points (reddish/light tones, circles) are 
for the RFIM. The right set of points (greenish/dark tones, squares) are for the DAFF. The large circles for the RFIM (large 
squares for the DAFF) represent our estimates of g{(,L/L) at T — Tc- The inset zooms into the important region (large box), 
where the Cardy Mapping should apply. Clearly, both data sets are different, suggesting that the RFIM and the DAFF do not 
share the same universality class. 



ior of the specific heat (critical exponent a). However, 
because the exponent is close to zero for both models, 
large system sizes are needed; sizes that arc currently 
not accessible via simulations. 

We conclude by cautioning researchers when using the 
equivalence of both models, in particular in experimental 
studies. 
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